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ABSTRACT 

This article presents a comprehensive framework for valuing financial instruments subject to credit 
risk. In particular, we focus on the impact of default dependence on asset pricing, as correlated default risk 
is one of the most pervasive threats in financial markets. We analyze how swap rates are affected by bilateral 
counterparty credit risk, and how CDS spreads depend on the trilateral credit risk of the buyer, seller, and 
reference entity in a contract. Moreover, we study the effect of collateralization on valuation, since the 
majority of OTC derivatives are collateralized. The model shows that a fully collateralized swap is risk¬ 
free, whereas a fully collateralized CDS is not equivalent to a risk-free one. 
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1 Introduction 


A broad range of financial instruments bear credit risk. Credit risk may be unilateral, bilateral, or 
multilateral. Some instruments such as, loans, bonds, etc, by nature contain only unilateral credit risk 
because only the default risk of one party appears to be relevant, whereas some other instruments, such as, 
over the counter (OTC) derivatives, securities financing transactions (SFT), and credit derivatives, bear 
bilateral or multilateral credit risk because two or more parties are susceptible to default risk. This paper 
mainly discusses bilateral and multilateral credit risk modeling, with a particular focus on default 
dependency, as correlated credit risk is one of the greatest threats to global financial markets. 

There are two primary types of models that attempt to describe default processes in the literature: 
structural models and reduced-form (or intensity) models. Many practitioners in the credit trading arena 
have tended to gravitate toward the reduced-from models given their mathematical tractability. They can 
be made consistent with the risk-neutral probabilities of default backed out from corporate bond prices or 
credit default swap (CDS ) spreads/premia. 

Central to the reduced-form models is the assumption that multiple defaults are independent 
conditional on the state of the economy. In reality, however, the default of one party might affect the default 
probabilities of other parties. Collin-Dufresne et al. (2003) and Zhang and Jorion (2007) find that a major 
credit event at one firm is associated with significant increases in the credit spreads of other firms. Giesecke 
(2004), Das et al. (2006), and Lando and Nielsen (2010) find that a defaulting firm can weaken the firms in 
its network of business links. These findings have important implications for the management of credit risk 
portfolios, where default relationships need to be explicitly modeled. 

The main drawback of the conditionally independent assumption or the reduced-form models is 
that the range of default correlations that can be achieved is typically too low when compared with empirical 
default correlations (see Das et al. (2007)). The responses to correct this weakness can be generally 
classified into two categories: endogenous default relationship approaches and exogenous default 
relationship approaches. 
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The endogenous approaches include the contagion (or infectious) models and frailty models. The 
frailty models (see Duffie et al. (2009), Koopman et al. (2011), etc) describe default clustering based on 
some unobservable explanatory variables. In variations of contagion or infectious type models (see Davis 
and Lo (2001), Jarrow and Yu (2001), etc.), the assumption of conditional independence is relaxed and 
default intensities are made to depend on default events of other entities. Contagion and frailty models fill 
an important gap but at the cost of analytic tractability. They can be especially difficult to implement for 
large portfolios. 

The exogenous approaches (see Li (2000), Laurent and Gregory (2005), Hull and White (2004), 
Brigo et al. (2011), etc) attempt to link marginal default probability distributions to the joint default 
probability distribution through some external functions. Due to their simplicity in use, the exogenous 
approaches become very popular in practice. 

Collateralization is one of the most important and widespread credit risk mitigation techniques used 
in derivatives transactions. According the ISDA (2012), 71% of all OTC derivatives transactions are subject 
to collateral agreements. The use of collateral in the financial markets has increased sharply over the past 
decade, yet the research on collateralized valuation is relatively sparse. Previous studies seem to turn away 
from direct and detailed modeling of collateralization (see Fuijii and Takahahsi (2012)). For example, 
Johannes and Sundaresan (2007), and Fuijii and Takahahsi (2012) characterize collateralization via a cost- 
of-collateral instantaneous rate (or stochastic dividend or convenience yield). Piterbarg (2010) regards 
collateral as a regular asset in a portfolio and uses the replication approach to price collateralized contracts. 

This paper presents a new framework for valuing defaultable financial instruments with or without 
collateral arrangements. The framework characterizes default dependencies exogenously, and models 
collateral processes directly based on the fundamental principals of collateral agreements. Some well- 
known risky valuation models in the markets, e.g., the CDS model, the risky interest rate swap (IRS) model 
(Duffie and Huang (1996)), can be viewed as special cases of this framework, when the default 
dependencies are ignored. 
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IRSs and CDSs arc two of the largest segments of the OTC derivatives market, collectively 
accounting for around two-thirds of both the notional amount and market value of all outstanding 
derivatives. Given this framework, we are able to analyze the value of IRSs with bilateral credit risk and 
look at how swap rates are affected by correlated default risk. Our study shows that counterparty default 
correlations have a relatively small impact on swap rates. Furthermore, we find that the value of a fully 
collateralized IRS is equal to the risk-free value. This conclusion is consistent with the current market 
practice in which market participants commonly assume fully collateralized swaps are risk-free. 

We also study the value of CDS contracts with trilateral credit risk and assess how spreads depend 
on the risk of the buyer, seller, and reference entity in a CDS contract. In general, a CDS contract is used 
to transfer the credit risk of a reference entity from one party to another. The risk circularity that transfers 
one type of risk (reference credit risk) into another (counterparty credit risk) within the CDS market is a 
concern for financial stability. Some people claim that the CDS market has increased financial contagion 
or even propose an outright ban on these instruments. 

The standard CDS pricing model in the market assumes that there is no counterparty risk. Although 
this oversimplified model may be accepted in normal market conditions, its reliability in times of distress 
has recently been questioned. In fact, counterparty risk has become one of the most dangerous threats to 
the CDS market. For some time now it has been realized that, in order to value a CDS properly, counterparty 
effects have to be taken into account (see ECB (2009)). 

We bring the concept of comvariance into the area of credit risk modeling to capture the statistical 
relationship among three or more random variables. Comvariance was first introduced to economics by 
Deardorff (1982), who used this measurement to correlate three factors in international trading. Furthermore, 
we define a new statistics, comrelation, as a scaled version of comvariance. Accounting for default 
correlations and comrelations becomes important in determining CDS premia, especially during the credit 
crisis. Our analysis shows that the effect of default dependencies on CDS premia from large to small is i) 
the default correlation between the protection seller and the reference entity, ii) the default comrelation, iii) 
the default correlation between the protection buyer and the reference entity, and iv) the default correlation 
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between the protection buyer and the protection seller. In particular, we find that the default 
comvariancc/comrelation has substantial effects on the asset pricing and risk management, which have 
never been documented. 

There is a significant increase in the use of collateral for CDS after the recent financial crises. Many 
people believe that, if a CDS is fully collateralized, there is no risk of failure to pay. Collateral posting 
regimes are originally designed and utilized for bilateral risk products, e.g., IRS, but there are many reasons 
to be concerned about the success of collateral posting in offsetting the risk of CDS contracts. First, the 
value of CDS contracts tends to move very suddenly with big jumps, whereas the price movement of IRS 
contracts is far smoother and less volatile than CDS. Second, CDS spreads can widen very rapidly. Third, 
CDS contracts have many more risk factors than IRS contracts. In fact, our model shows that full 
collateralization cannot eliminate counterparty risk completely for a CDS. 

This article also shows that the pricing process of a defaultable instrument normally has a backward 
recursive nature if the payoff can be positive or negative. Accordingly, we propose a backward induction 
approach for risky valuation. In contrast to the popular recursive integral solution (see Duffie and Fluang 
(1996)), our backward induction method significantly simplifies the implementation. One can make use of 
the well-established algorithms, such as lattice/tree and regression-based Monte Carlo, to price a defaultable 
instrument. 

The rest of this paper is organized as follows: Pricing bilateral defaultable instruments is elaborated 
on in Section 2; valuing multilateral defaultable instruments is discussed in Section 3; the conclusions are 
presented in Section 4. All proofs and some detailed derivations are contained in the appendices. 

2 Pricing Financial Instruments Subject to Bilateral Credit Risk 

We consider a filtered probability space (Q , J 7 , \p t [> 0 , P ) satisfying the usual conditions, 
where Q denotes a sample space, J 7 denotes a a -algebra, P denotes a probability measure, and 
ULo denotes a filtration. 
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In the reduced-form approach, the stopping (or default) time r, of firm i is modeled as a Cox arrival 
process (also known as a doubly stochastic Poisson process) whose first jump occurs at default and is 
defined by, 

r, =inf {t-\ o hi (s, Z s )ds>H t j (1) 

where h t (t) or h i {t,Z l ) denotes the stochastic hazard rate or arrival intensity dependent on an exogenous 
common state Z t , and //, is a unit exponential random variable independent of Z,. 

Dependence between the default times is only introduced by the dependence of the intensity /r, (/) 
on a common process Z,. Consequently, conditional on the path of Z f , defaults arc independent, which is 

the reason why this setup is also often called the conditional independence setup. 

It is well-known that the survival probability from time t to s in this framework is defined by 

(t,s):=P j (r> s\z>t,Z t ) = exp^- J h i (u)du) (2a) 

The default probability for the period (t, 5 ) in this framework is given by 

:= Pfiz < s I t > t,Z t ) = l-p ; (f,s) = l-expf-| h^^duj (2b) 

Three different recovery models exist in the literature. The default payoff is either i) a fraction of 
par (Madan and Unal (1998)), ii) a fraction of an equivalent default-free bond (Jarrow and Turnbull (1995)), 
or iii) a fraction of market value (Duffie and Singleton (1999)). The whole course of the recovery 
proceedings under the Bankruptcy and Insolvency act is a complex process that typically involves extensive 
negotiation and litigation. No model can fully capture all aspects of this process so, in practice, all models 
involve trade-offs between different perspectives and views. In general, the choice for a certain recovery 
assumption is based on the legal structure of an instrument to be priced. For example, the recovery of market 
value (RMV) assumption is well matched to the legal structure of an IRS contract where, upon default 
close-out, valuation will in many circumstances reflect the replacement cost of the transaction, whereas the 
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best default recovery assumption for a CDS is that the claim made in the event of the reference default 
equals a fraction of the face value of the underlying bond. 

There is ample evidence that corporate defaults are correlated. The default of a firm’s counterparty 
might affect its own default probability. Thus, default correlation/dependence arises due to the counterparty 
relations. 

Two counterparties are denoted as A and B. The binomial default rule considers only two possible 
states: default or survival. Therefore, the default indicator Y j for party j (j-A, B ) follows a Bernoulli 

distribution, which takes value 1 with default probability q jt and value 0 with survival probability p , i.e., 
P{Yj = 0} = pj and P{ Y s = 1} = q . The marginal default distributions can be determined by the reduced- 

form models. The joint distributions of a multivariate Bernoulli variable can be easily obtained via the 
marginal distributions by introducing extra correlations. 

Consider a pair of random variables (Y A ,Y B ) that has a bivariate Bernoulli distribution. The joint 


probability representations are given by 

Poo ■= P(Y a =0,Y b = 0) = p A p B + cr^ (3a) 

Po\ ■= p ( y A = °’ Y b= 0 = PaC/b - °AR (3b) 

P\o : = p ( y A = b y B = 0) = 4 a Pk ~ °ak (3c) 

Pi 1 : = P ( y A =IY B =\) = q A P B + °AB (3d) 


where E{Y J ) = q j , a) = p i q t , and <j ab := e[(Y a - q A )(Y B - q B )\ = p AB o A a B = p AB Jq A P A q B P B where p AB 
denotes the default correlation coefficient, and a AB denotes the default covariance. 

A critical ingredient of the pricing of a bilateral defaultable instrument is the default settlement 
rules. There are two rules in the market. The one-way payment rule was specified by the early International 
Swap Dealers Association (ISDA) master agreement. The non-defaulting party is not obligated to 
compensate the defaulting party if the remaining market value of the instrument is positive for the defaulting 
party. The tw’o-way payment rule is based on current ISDA documentation. The non-defaulting party will 
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pay the full market value of the instrument to the defaulting party if the contract has positive value to the 
defaulting party. 

1.1 Risky valuation without collateralization 

Consider a defaultable instrument that promises to pay a X T from party B to party A at maturity 
date T, and nothing before date T. The payoff X T may be positive or negative, i.e. the instrument may be 
either an asset or a liability to each party. All calculations are from the perspective of party A. 

We divide the time period (t , T) into n very small time intervals (At) and use the approximation 
exp(y)~ 1 + y provided that y is very small. The survival and the default probabilities for the period (t, 
t + At ) are given by 

p(t) := p(t, t + At) = exp (- h(t)At) «1 - h(t)At (4a) 

q(t) := q(t , t + At) = 1 - exp (- h(t)At) « h(t)At (4b) 

Suppose that the value of the instrument at time t + At is V(t + At) that can be an asset or a liability. 

There are a total of four ( 2 2 = 4) possible states shown in Table 1. 

The risky value of the instrument at time t is the discounted expectation of all the payoffs and is 
given by 

V(t) = £{exp(- r(t)At)[l ni+At) > 0 {p 00 (t) + <p B (t)p 0l (t)+ p B (t)p l0 (t) + <p AB (t)p n (t))V(t + At)\J t 

ly(/+AO<0 (Poo (0 + (Pa iOPoi (0+ (P A it)p l0 ( t) + (p AB (t)p l , (t) )V (t + At)\j t ]} ( 5a ) 

« £{exp[- {r(t) + 

V(f+Af)>0 l B (t) + ^-V(t+At)<0^A (t))\r]k ( t + At)\f t } = .E{exp(- g(t)At)V ( t + At)\f t } 

where 

g(t) = r(t)+ 

^-V(t+At)>0^B (0 (t+At)<0^A (t) (5b) 

l b^) = (l -cp B (t))h B (t) + (l ~(p B (t))h A (t) - (l ~<p B (t) - <p B (t) + (,o AB ( t))p AB (t)yjh A ( t)h B ( t) (5c) 

l a i f ) = (l-<A it))h A i0 + (l ~9 a it)) h B it) - (l ~(Pa i f ) ~ (P A i f ) + <Pab it))p AB it)ylh A it)h B it) (5d) 

where \ Y is an indicator function that is equal to one if Y is true and zero otherwise, £ !• \PF, j is the 
expectation conditional on the f t , r(t) is the risk-free short rate, and tp t is the recovery rate. 
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The pricing equation above keeps terms of order At. All higher order terms of At are omitted. 
Similarly, we have 

V(t + At) = £jexp(- g(t + At)At)V(t + 2At)\jT, +Al } (6) 

Note that exp(-g(t)Af) is J/ +A/ -measurable. By definition, an 5Fu\t -measurable random 

variable is a random variable whose value is known at time t + At. Based on the taking out what is known 
and tower properties of conditional expectation, we have 

V(t) = £{exp(- g(t)At)V(t + At)\jF,}= £{exp(-g(t)At)f?[exp(- g(t + At)At)V(t+ 2At)\jF t+At ] 

= ^{exp(- X ) =0 g(! + iAt)At)jv (t + 2At)\f t } 

By recursively deriving from t forward over T where V(T) = X t and taking the limit as At 
approaches zero, we obtain 


V(t) = E{G(f,T)X T \jF t }=E{ exp|~-j^ g(u)du 


At f. 


(8) 


We may think of G(t,T) as the bilateral risk-adjusted discount factor and g(u) as the bilateral risk- 
adjusted short rate. Equation (8) has a general form that applies in a particular situation where we assume 
that parties A and B have independent default risks, i.e. p AB = 0 and <p AB = 0. Thus, we have: 


V(t) = E{G(t, T)X t \JF, }= Ejexpf g(u)du 


At AFt 


(9a) 


where 


g(u) = r(u) + l v{u) > 0 l B (u) + l vw<0 l A (u) 

(9b) 

l B (u) = (l -cp B ( u))h B (u) + i)--(p B 0 uj)h A (u) 

(9c) 

l A (i u ) = (l - (p A ( u))h A (u) + (l - (p A ( u))h B (i u ) 

(9d) 


Equation (9) is the same as equation (2.5’) in Duffie and Huang (1996). 

In theory, a default may happen at any time, i.e., a risky contract is continuously defaultable. This 
Continuous Time Risky Valuation Model is accurate but sometimes complex and expensive. For simplicity. 
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people sometimes prefer the Discrete Time Risky Valuation Model that assumes that a default may only 
happen at some discrete times. A natural selection is to assume that a default may occur only on the payment 
dates. Fortunately, the level of accuracy for this discrete approximation is well inside the typical bid-ask 
spread for most applications (see O’Kane and Turnbull (2003)). From now on, we will focus on the discrete 
setting only, but many of the points we make arc equally applicable to the continuous setting. 

If we assume that a default may occur only on the payment date, the risky value of the instrument 


in a discrete-time setting is given by 

V(t) = E{D(t,T)[l XT > 0 { Poo (t,T) + (p B (T)p ol (t,T)+p B (T)p w (t,T) + <p AB (T)p n (t,T))X T \f t 

+ix r <o(Poo^r)+^(r)Po 1 (/,r)+^(r)p 10 (fT)+^(r)p 1I (^r))x r |jJ n0a ) 

= e[d(1, T)( i X t> _ 0 k B ( t,T) +1 Xr<0 k A 0 1, T))x t \f t ] = e{k(!,T)X t \? t ) 

where 

k B (t,T) = p B (t,T)p A (t,T) + <p B (T)q B (t,T)p A (t,T) + <p B (T)p B (t,T)q A (t,T) 

+ <Pab (F t )<1a (F r ) + ° AB (?’ T){l- l p B (T)-ip B (J’) + ip AB (T)) Ui 11 


k A (t,T) = p B {t,T)p A (t,T) + (p A ( T)q A (f, T)p B (t,T) + <p A ( T)p A ( t , T)q B (f, T) 

+ <Pab b (0 t )<1a (0 T ) + ° ab (0 T ) 0 - <Pa (T) ~ <Pa ( t ) + <Pab i T )) 

where D(t,r) denotes the stochastic risk-free discount factor at t for the maturity T given by 


D(t,T) = exp 



(10c) 


(10d) 


We may think of K(t,T) as the risk-adjusted discount factor, and k A (t,T ) and k B (t,T ) as the 
adjustment factors. Equation (10) tells us that the bilateral risky price of a single-payment instrument can 
be expressed as the present value of the payoff discounted by a risk-adjusted discount factor that has a 
switching-type dependence on the sign of the payoff. 

Equation (10) can be easily extended from one-period to multiple-periods. Suppose that a 
defaultable instrument has m cash flows. Let the m cash flows be represented as X m with payment 

dates 7j,..., T m . Each cash flow may be positive or negative. We have the following proposition. 
Proposition 1: The risky value of the multiple-payment instrument is given by 
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(Ha) 


v «-ir.i^(n'‘W r i.wkwl 

where t = T 0 and 

K ( T j’ T j +ll = D(Tj , Tj+I ) (l( x j+i +V(T J+l ))>0^B (Tj ’ Tj+i ) + 1(X J+1 +V (7) +1 »<o ^ a (Tj » Tj+i ) j (1 lb) 

where k A (Tj, T n i ) and Ji^(Tj ,7j^j) ate defined in Ec/uation (10). 

Proof: See the Appendix. 

From Proposition 1, we can see that the intermediate values are vital to determine the final price. 
For a payment interval, the current risky value has a dependence on the future risky value. Only on the final 
payment date T m , the value of the instrument and the maximum amount of information needed to determine 


the risk-adjusted discount factor are revealed. This type of problem can be best solved by working 
backwards in time, with the later risky value feeding into the earlier ones, so that the process builds on itself 
in a recursive fashion, which is referred to as backward induction. The most popular backward induction 
valuation algorithms are lattice/tree and regression-based Monte Carlo. 

1.2 Risky valuation with collateralization 

Collateralization is the most important and widely used technique in practice to mitigate credit risk. 
The posting of collateral is regulated by the Credit Support Annex (CSA) that specifies a variety of terms 
including the threshold, the independent amount, and the minimum transfer amount (MTA), etc. The 
threshold is the unsecured credit exposure that a party is willing to bear. The minimum transfer amount is 
the smallest amount of collateral that can be transferred. The independent amount plays the same role as 
the initial margin (or haircuts). 

In a typical collateral procedure, a financial instrument is periodically marked-to-market and the 
collateral is adjusted to reflect changes in value. The collateral is called as soon as the mark-to-market 
(MTM) value rises above the given collateral threshold, or more precisely, above the threshold amount plus 
the minimum transfer amount. Thus, the collateral amount posted at time t is given by 


C(t) = 


(V(t)-H(t) 

[0 


if V(t) > H(t) 
otherwise 


( 12 ) 
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where H{t) is the collateral threshold. In particular, H(t) = 0 corresponds to full-collateralization 2 ; // > 0 
represents partial/under-collateralization; and H< 0 is associated with over-collateralization. Full 
collateralization becomes increasingly popular at the transaction level. In this paper, we focus on full 
collateralization only, i.e., C(t) =V(t). 

The main role of collateral should be viewed as an improved recovery in the event of a counterparty 
default. According to Bankruptcy law, if there has been no default, the collateral is returned to the collateral 
giver by the collateral taker. If a default occurs, the collateral taker possesses the collateral. In other words, 
collateral does not affect the survival payment; instead, it takes effect on the default payment only. 

For a discrete one-period (t, u) economy, the posted collateral at time t is C(t) =V(t). At time u, 
there are several possible states: i) Both A and B survive. The instrument value is equal to the market value 
V(u)\ ii) Either or both counterparties A and B default. The instrument value is the future value of the 
collateral, i.e., C(u) = V(t)l D(tji) where we consider the time value of money only. Since the majority of 
the collateral is cash according to ISDA (2012), it is reasonable to consider the time value of money only 
for collateral assets. The large use of cash means that collateral is both liquid and not subject to large 
fluctuations in value. The value of the collateralized instrument at time t is the discounted expectation of 
all the payoffs and is given by 

V(t ) - E[D(t , u)(p 0Q (t, u)V (w) + p 0l (t, m)C(w) + p l0 (t, m)C(m) + p u (t,u)C{u))\f t ] 

= e[d(T u) p 00 (t,u)V (, u ) | ft ]+ [l - E(p 00 (t, u) \j t )]v (0 

or 

V(t) = E[D(t, s)p 00 (t, s)V (s) | J'f ]/ e(p 00 (t, s) \jf t ) (13b) 


2 There are three types of collateralization: Full-collateralization is a process where the posting of collateral 
is equal to the current MTM value. Partial/under-collateralization is a process where the posting of collateral 
is less than the current MTM value. Over-collateralization is a process where the posting of collateral is 
greater than the current MTM value. 
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If we assume that default probabilities arc uncorrelated with interest rates and payoffs 3 , we have 

V{t) = E[D{t,s)V{s)\y r t \ (14) 

Equation (14) is the formula for the risk-free valuation. Thus, we have the following proposition. 
Proposition 2: If a bilateral risky instrument is fully collateralized, the risky value of the instrument is equal 
to the risk-free value, as shown in equation (14). 

Since an IRS is a typical bilateral risky contract. Proposition 2 squares with the results of Johannes 
and Sundaresan (2007), and is also consistent with the current market practice in which market participants 
commonly assume fully collateralized swaps are risk-free and it is common to build models of swap rates 
assuming that swaps are free of counterparty risk. 

1.3 Numerical results 

To study the impact of bilateral credit risk, we have selected a new 10-year fixed-for-floating IRS 
with a quarterly payment frequency. Two counterparties are denoted as A and B. Counterparty A pays a 
fixed rate and counterparty B pays a floating rate. All calculations are from the perspective of party A. 

The current (spot) market data are shown in Table 2, provided by FinPricing (2019). At the time 
the contract is entered into, there is no advantage to either party. By definition, a swap rate is a fixed rate 
that makes the market value of a given swap at initiation zero. The risk-free swap rate can be easily 
calculated as 0.03433. 

Since the payoff of an IRS is determined by interest rates, we need to model the evolution of the 
floating rates. Interest rate models are based on evolving either short rates, instantaneous forward rates, or 
market forward rates (e.g., the LIBOR Market Model (LMM)). Since both short rates and instantaneous 
forward rates are not directly observable in the market, the models based on these rates have difficulties in 
expressing market views and quotes in term of model parameters, and lack agreement with market valuation 


3 Moody’s Investor’s Service (2000) presents statistics that suggest that the correlations between interest 
rates, default probabilities, and recovery rates are very small and provides a reasonable comfort level for 
the uncorrelated assumption. 
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formulas for basic derivatives. On the other hand, the object modeled under the LMM is market-observable. 
It is also consistent with the market standard approach for pricing caps/floors using Black’s formula. They 
are generally considered to have more desirable theoretical calibration properties than short rate or 
instantaneous forward rate models. Therefore, we choose the LMM lattice proposed by Xiao (2011) for 
pricing the defaultable IRS. We also implement the Hull-White trinomial tree to verify the results and 
ensure robustness of the valuation. This paper, however, only reports the results produced by the LMM 
lattice. 

We first assume that i) counterparties A and B have independent default risks; ii) the hazard rates 
are deterministic; and iii) both parties have a constant recovery of 60%. We use the LMM to evolve the 
interest rates and then price the risky IRS according to Proposition 1. The risky swap rates are computed 
and shown in Table 3. 

From Table 3, we derive the following conclusions: First, a fixed-rate payer with lower credit 
quality (higher credit risk) pays a higher fixed rate. Second, a credit spread of about 100 basis points 
translates into a swap spread of about 1.3 basis points. Finally, the credit impact on swap rates is 
approximately linear within the range of normally encountered credit quality. This confirms the findings of 
Duffie and Huang (1996). Intuitively, a risk-free floating-rate payer demands a higher fixed rate if the fixed- 
rate payer has a lower credit score. 

We next present some new results. Assume that party A has an ‘A+300bps’ credit quality, i.e., a 
‘300 basis points’ parallel shift in the A-rated credit spreads, and party B has an ‘A’ credit quality. The 
risky swap rate with asymmetric credit qualities is calculated as 0.03436. 

Assume <p AB = 0.5. The effect of the default correlation p AB on swap rate is shown in Figure 1. We 
can draw the following conclusions from the results: First, the counterparty default correlation and the swap 
rates have a negative relation, i.e., a negative sensitivity of swap rates to changes in counterparty default 
correlation is obtained. Second, the graph suggests an almost linear relationship between the swap rates and 
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the default correlation. Finally, the impact of the default correlation is modest (e.g., in the range of [-2, 2] 
basis points). 

Bilateral credit risk modeling is probably the simplest example involving default dependency, but 
it shows several essential features for modeling correlated credit risk, which will help the reader better 
understand the increasingly complex cases in the following section. 

2 Pricing Financial Instruments Subject to Multilateral Credit Risk 

The interest in the financial industry for the modeling and pricing of multilateral defaultable 
instruments arises mainly in two respects: in the management of credit risk at a portfolio level and in the 
valuation of credit derivatives. Central to the pricing and risk management of credit derivatives and credit 
risk portfolios is the issue of default relationships. 

Let us discuss the three-party case first. A CDS is a good example of a trilateral defaultable 
instrument where the three parties are counterparties A, B and reference entity C. In a standard CDS contract 
one party purchases credit protection from another party, to cover the loss of the face value of a reference 
entity following a credit event. The protection buyer makes periodic payments to the seller until the maturity 
date or until a credit event occurs. A credit event usually requires a final accrual payment by the buyer and 
a loss protection payment by the protection seller. The protection payment is equal to the difference between 
par and the price of the cheapest to deliver (CTD) asset of the reference entity on the face value of the 
protection. 

A CDS is normally used to transfer the credit risk of a reference entity between two counterparties. 
The contract reduces the credit risk of the reference entity but gives rise to another form of risk: 
counterparty risk. Since the dealers arc highly concentrated within a small group, any of them may be too 
big to fail. The interconnected nature, with dealers being tied to each other through chains of OTC 
derivatives, results in increased contagion risk. Due to its concentration and interconnectedness, the CDS 
market seems to pose a systemic risk to financial market stability. In fact, the CDS is blamed for playing a 
pivotal role in the collapse of Lehman Brothers and the disintegration of AIG. 
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For years, a widespread practice in the market has been to mark CDS to market without taking the 
counterparty risk into account. The realization that even the most prestigious investment banks could go 
bankrupt has shattered the foundation of the practice. It is wiser to face frankly the real complexities of 
pricing a CDS than to indulge in simplifications that have proved treacherous. For some time now it has 
been realized that, in order to value a CDS properly, counterparty effects have to be taken into account. 

The default indicator K ; for firm j (j = A or B or C) follows a Bernoulli distribution, which takes 

value 1 with default probability q } , and value 0 with survival probability p . The joint probability 


representations of a trivariate Bernoulli distribution (see Teugels (1990)) are given by 

Pooo := PO' A =0 ,Y b =0 ,Y C =0) = p A p B Pc + Pc u ab + Pb (7 ac + Pa (T bc ~®abc (15a) 

Fioo :=y> ^ y a = 1 ~ y b = 0,T C =0) = <1aPbPc ~ Pc a AB ~ Pb^ac +( 1a <j bc + @abc (15b) 

Foio :=F(T a =0 ,Y b =1,F C =0) = p A q B Pc ~ Pc^ab +c Lb°ac ~Pa^bc +0 abc (15c) 

FoOl := ^ Y A = 0 ’ Y R = ®’ Y C = 1 PaPbQc +C lc <7 AB ~ PB a AC ~ PA a BC + @abc (15d) 

Puo := ^ y a = 1 > y b = l’lc = H) = QaQbPc + Pc° ab^Qb 17 ac ~Pa°bc ~®abc (15e) 

Fioi -—Pi}' A = I- Y b =0 ,Y C =1 ) = QaPbQc ~ c 1c°ab + Pb <j ac ~ c 1a (7 bc ~®abc (15f) 

Foil == p ( y a =^ y b =1 > Y c =1 ) = Pa<Ib<Ic -Qc^ab ~<1b°ac + Pa°bc ~ 9 abc (15g) 

Fm :=y> ( y a ~ yy b =yy c = 1) =( Ia c 1b c 1c +< 1c <j ab +c 1b (7 ac + c 1a (T bc + @abc (15h) 

where 

®abc := e{( y a -<1a)( y b ~ c Ib)( y c ~ c lc)) (15i) 


Equation (15) tells us that the joint probability distribution of three defaultable parties depends not 
only on the bivariate statistical relationships of all pair-wise combinations (e.g., <r..) but also on the 

bivariate statistical relationship (e.g., 0 ABC ). 0 ABC was first defined by Deardorff (1982) as comvariance, 

who use it to correlate three random variables that are the value of commodity net imports/exports, factor 
intensity, and factor abundance in international trading. 
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We introduce the concept of comvariance into credit risk modeling arena to exploit any statistical 
relationship among multiple random variables. Furthermore, we define a new statistic, comrelation, as a 
scaled version of comvariance (just like correlation is a scaled version of covariance) as follows: 
Definition 1: For three random variables X A , X B ,and X c ,let ju A , p B , and /i c denote the means of X A 
, X B ,and X c . The comrelation of X A , X B ,and X c is defined by 


C ABC - 


e[(X a - M a)(X b -M b )(X c -v c )\ 
3 fE\ X A xE\X B - ju b | 3 xE\X c -ju c 


(16) 


According to the Flolder inequality, we have 

\e{{X a -H a )(X b -m b )(X c -fi c ))\< E\(X a -ju a )(X b -ju b )(X c -V c )\<^E\X a - Ma \ 3 xE\X b - Mb \ 3 x^X c - Mc \ 3 (17) 

Obviously, the comrelation is in the range of [-1, 1], Given the comrelation, Equation (15i) can be 
rewritten as 


e ABC := E((Y A - q A )(Y B - q B )(X c - q c j) = Cabc^E\Y a - q A f x E\Y b - q B f x E\Y C - q c f 

I -------— (18) 

= £ABCy^\P A c 1a(P A + c Ia)Pb c Ib^Pb + C 1B^ P c c lc(p c +c lc\ 

where E{Y j ) = q. and E\Y j -q.| 3 = py/fp) I q]), j=A, B, or C. 


If we have a series of n measurements of X A , X B , and X c written as X Ai , X Bi and x Ci where / - 
1 ,2 ,the sample comrelation coefficient can be obtained as: 


CABC - ' 


X"= 1 ( x At - Ba )( x Bi ~ Be )(*ci ~ Be ) 


(19) 


-i / v - ' n | 1 3 ' n | 1 3 ' n i 1 3 

V2w= 1 1 x Ai - Ba | x Li= 1 \ x Bi -Be | x X, i=l \x a - Be | 

More generally, we define the comrelation in the context of n random variables as 
Definition 2: For n random variables X lt X 2 ,..., X n , let denote the mean of X ; where i=l,..,n. The 
comrelation of X ,, X 2 , X n is defined as 

E[(X 1 -ju l )(X 2 -B 2 )-(X n - Mn )] 


C\2...n ~ 


E\Xi-Bi\ n xE\X 2 -Bi\ n -xE\X n -ju n \ n 


( 20 ) 
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The correlation is just a specific case of the comrelation where n - 2. Again, the comrelation Cl2...n 


is in the range of [-1, 1] according to the Holder inequality. 

2.1 Risky valuation without collateralization 

Recovery assumptions arc important for pricing credit derivatives. If the reference entity under a 
CDS contract defaults, the best assumption, as pointed out by J. P. Morgan (1999), is that the recovered 
value equals the recovery rate times the face value plus accrued interest 4 . In other words, the recovery of 
par value assumption is a better fit upon the default of the reference entity, whereas the recovery of market 
value assumption is a more suitable choice in the event of a counterparty default. 

Let valuation date be t. Suppose that a CDS has m scheduled payments. Let each payment be 
represented as X t =-sNS(T j _ l ,T i ) with payment dates T l ,...,T m where i=l,„,m, 8{T j _ l ,T j ) denotes the 
accrual factor for period (T j ,, 7 ), N denotes the notional/principal, and s denotes the CDS premium. Party 

A pays the premium/fee to party B if reference entity C does not default. In return, party B agrees to pay the 
protection amount to party A if reference entity C defaults before the maturity. We have the following 
proposition. 

Proposition 3: The value of the multiple-paymen t CDS is given by 


vw=s:i,d(n;>v J .r-)k 1 |j,j + z:: 1 d(n;>(r.r., 


jkl!7: l TjRlT_ r Tj\T 1 \ (21a) 


where t = T 0 and 

0(Tj,T j+l ) = l(v (r . + j +x . +l ^o^s(^>^u) + l(y(r J+ n+x /+1 ko^A(^,7)+i) (21b) 

h (Tj , T j+1 ) = { Pa (ITj , T J+l )p B (Tj , T j+1 ) Pc (T j , T j+1 ) + q A (7\, T j+1 )p B {T j , T J+1 )p c (T j , T j+1 )cp A (T ]+1 ) 

+ P A ( Tj , T j+1 )q B (T p T j+] ) Pc (T., T j+l )<p A (T '. +1 ) + q A (T j , T '. +1 )q B (T ] ,T J+l )p c (T j , T j+1 )(p AB (T j+1 ) 

+ Pc (Tj, T j+1 )a AB C Tj , T J+1 )(l -<p A (T j+1 ) - <p A (T j+1 ) + cp AB (T j+1 )) (21c) 

+ VAC (T i . T j + i)[pb (Tj . T i +l )(> - (Pa (Tj + 1 ))+ (1b (Tj . T i+ i){<P a ( T j+ i ) “ Vab (T ) + i ))] 

+ cr B c ( T J . T )+i )[pa (Ti . T .i +l )(' - (Pa (T) + i ))+ (,l a (Tj . T j+ i ){<Pa ( T j+ i ) “ <Pab C T J+1 ))] 

+ Oabc (Tj , T j+l ){-1 + ip A (T j+I ) - cp AB (T j+1 ) + cp A (T j+1 ))}d{T ] , T ]+i ) 


4 In the market, there is an average accrual premium assumption, i.e., the average accrued premium is half 
the full premium due to be paid at the end of the premium. 
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K (Tj , T j+l )={p A aTj , T j+l )p B (T,,T j+1 ) p c (T j , T j+1 ) + q A {T p T J+l )p B (T j , T j+1 )p c (T j ,7\ +1 )cp B (T J+I ) 

+ P a Wj > Tj+i )q B (Tj, Tj+i )p c (Tj, Tj +] )cp B (Tj+i) + q A (Tj, T j+l )q B (Tj, T j+l )p c (Tj, T J+1 )(p AB (Tj +l ) 

+ Pc (Tj , Tj+i )a AB (Tj , Tj+i )(l -<p B (Tj+i ) - cp B (T j+1 ) + cp AB (T j+1 )) (2 Id) 

+ o’ac (Tj,Tj+i)[p B (Tj , T j+l )(l — q> B (Tj-+i ))+ q B (Tj,Tj +l )[(p B (7j +1 )~ Vab (^j+i ))] 

+ a BC (Tj , T j+1 )[p A (Tj ,T J+ i )(l - cp B (Tj )) + q A (7\, T j+1 )(<p B (T j+l ) - (p AB (T j+1 ))] 

+ 6(Tj , Tj+i )(-1 + cp B (Tj+i ) - cp AB (Tj+i ) + cp B (Tj+i ))}D(Tj , T j+1 ) 

TKTj , Tj+i ) = [p A (Tj , Tj+i ) p B (Tj , Tj+i )q c (Tj , 7j +1 ) + q A (Tj , T J+l )p B (Tj , 7j +1 ) 9c (7., T J+1 )<p B (7j +1 ) 

+ Pa (2)>2) +1 )q B (Tj,Tj+i )q c (Tj , 7j +1 )^> g (7j. +1 ) + q A (Tj,T j+1 )q B (Tj, 7\ +1 )r/ c (Tj,7) +1 )^ Ag (7) +1 ) 

+ q c (Tj,Tj+l )a AB (Tj ,r. +1 )(l - cp B (Tj+l) - <p B (T j+1 ) + ^ Ag (7j +1 )) (21 e) 

- ^ac . r, +1 )[p b (r,, r. +1 )(i - ^ g (r, +1 ))+ q B (Tj , r. +1 ){<p B (r. +1 ) - ^ Ag (r. +1 ))] 

— cx gc (Tj,Tj+i )\p A (Tj,Tj+i )(l — <p B (Tj+i))+ q A (T ) ,Tj+i)(<p B ( T J+1 ) — <p AB (7) +1 ))] 

0ABC (Tj > Tj + 1 )(l - (r, +I ) + <Pab (T y+1 ) - cp B (Ti+SfjDVj , r. +I ) 

where 7?(T.,T. +1 ) = (iv(l-^ c (r. +1 ))-«(r.,r. +1 )), a(T j ,T j +i) = sN8(T s ,T)/ 2, and X,. = -sNS<X j ,T J + i )- 


Proof: See the Appendix. 

We may think of 0(t,T) as the risk-adjusted discount factor for the premium and Tl(t,T) as the 
risk-adjusted discount factor for the default payment. Proposition 3 says that the pricing process of a 
multiple-payment instrument has a backward nature since there is no way of knowing which risk-adjusted 
discounting rate should be used without knowledge of the future value. Only on the maturity date, the value 
of an instrument and the decision strategy are clear. Therefore, the evaluation must be done in a backward 
fashion, working from the final payment date towards the present. This type of valuation process is referred 
to as backward induction. 

Proposition 3 provides a general form for pricing a CDS. Applying it to a particular situation in 
which we assume that counterparties A and B are default-free, i.e., p f = 1, qj = 0, p kl =0, and g ABC = 0 , 
where j=A or B and k, l-A, B, or C, we derive the following corollary. 

Corollary 1: If counterparties A and B are default-free, the value of the multiple-payment CDS is given by 


=z^[(n ;:o ^ 1-^ J + z z, ^[(n ;: 2 o ^^ >i^ 


( 22 ) 


where 0(7)^,7)) = DCf i,T)p c (T ,, 7j ); 0(7^ , 7) ) = D(T P1 , 7} )q c (T pl , 7) ). 
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The proof of this corollary becomes straightforward according to Proposition 3 by setting p kl =0, 
<p AB = 0, W = 0 , Pj = 1, cjj = o , PcitJi) = YX;l {) p(T K J K+ 0 , and = YX^ D{T ^ T ^ ■ 

If we further assume that the discount factor and the default probability of the reference entity are 
uncorrelated and the recovery rate cp c is constant, we have 

Corollary 2: Assume that i) counterparties A and B are default-free, ii) the discount factor and the default 
probability of the reference entity are uncorrelated; Hi) the recovery rate cp c is constant; the value of the 

multiple-payment CDS is given by 

where P(t,f ) = E\p(t,T i )\J\[ denotes the bond price, p c it,T i ) = E\p c {t,T j )\ff r t \, q c {t,T i ) =1- p c {t,T i ), 
pit, 7}_, )q(Ti-\ ,T t ) = pit, 7)_!) - pit, 7j ). 

This corollary is easily proved according to Corollary 1 by setting t[a') / | j^]= /i'[A'| j^]/r[) / | j^] when 

X and Y are uncorrelated. Corollary 2 is the formula for pricing CDS in the market. 

Our methodology can be extended to the cases where the number of parties n > 4. A generating 
function for the (probability) joint distribution (see details in Teugels (1990)) of n-variate Bernoulli can be 
expressed as 



where ® denotes the Kronecker product; p < " ) = {pf’ j and a-"' 1 = joy" 1 1 are vectors containing 2 


components 


: p[ n) = Pk„k 2 ,...,k n . k = l+ T!U k i 2 ‘ 1 ’ k i e {°’ 1 }; 


2.2 Risky valuation with collateralization 

According to the ISDA (2012), almost all CDSs are fully collateralized. Many people believe that 
full collateralization can eliminate counterparty risk completely for CDS. 

Collateral posting regimes are originally designed and utilized for bilateral risk products, e.g., IRS, 
but there are many reasons to be concerned about the success of collateral posting in offsetting the risks of 
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CDS contracts. First, the values of CDS contracts tend to move very suddenly with big jumps, whereas the 
price movements of IRS contracts arc far smoother and less volatile than CDS prices. Second, CDS spreads 
can widen very rapidly. The amount of collateral that one party is required to provide at short notice may, 
in some cases, be close to the notional amount of the CDS and may therefore exceed that party’s short-term 
liquidity capacity, thereby triggering a liquidity crisis. Third, CDS contracts have many more risk factors 
than IRS contracts. 

We assume that a CDS is fully collateralized, i.e., the posting of collateral is equal to the amount 
of the current MTM value: C(t) = V(t ) . For a discrete one-period (t, u) economy, there are several possible 

states at time u : i) A, B, and C survive with probability p 000 . The instrument value is equal to the market 
value V(u ); ii) A and B survive, but C defaults with probability p QQl . The instrument value is the default 
payment R(u) ; iii) For the remaining cases, either or both counterparties A and B default. The instrument 
value is the future value of the collateral V(t)/ D(t,u) (Here we consider the time value of money only). 
The value of the collateralized instrument at time t is the discounted expectation of all the payoffs and is 
given by 

V ( t ) = E{D(t, u)[p 000 (t, u)V ( u) + p 001 (t, u)R(u) 

+ ( p m (t , a) + p 0l0 (t, u) + p l 10 (t, u) + p m (t, u) + p 0l ,(i, u) + p l n (t, u)) V ( t ) / D(t, u )]\f t } (25a) 

= E {[ D(t , u)(p 000 (t, u)V (s) + p 0(n ( t , u)R(u))+ (l - p m) (t, u) - p 0Ql (t,u))V(t)]\f t } 

or 

(f’ u )Pb (T u ) + &AB (T 

= E^D(t,u)(p c (t, u)V (u) + q c (t, u)R(u)){p A ( t , u)p B (i, u ) + a AB (i, u )) (25b) 

+ D(t, u)(V (m) - R(u))(p B (t, u)a AC (i, ii) + p A ( t,u)a BC (i, u ) - d ABC (i,«))] f t } 

If we assume that (p A (t,u)p B (t,u) + a AB (t,u)) and D(t,u)(p c (t,u)V(u) + q c (t,u)R(u)) are 

uncorrelated, we have 

V(t) =V F (t) +% ABC (t,u)h// AB (t,u) (26a) 

where 
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V F (t) = E{D(t,u)[p c (t,u)V(u) + q c (t,u)R(u)\\jy} 


(26b) 


W.AB (C M ) — e\\ p A (t, u)p B (t, u) + O' AB (t,u)]\jT t } (26c) 

^ABC (0 M ) = [0(t, U) (p fj (t, U)<J AC (t, U) + p A (t,U.)c7 BC (t,U) — 6 abc (t,u)){V(u)-R(u)]]\jr t } (26d) 

The first term V f (?) in equation (26) is the counterparty-risk-free value of the CDS and the second 
term is the exposure left over under full collateralization, which can be substantial. 

Proposition 4: If a CDS is fully collateralized, the risky value of the CDS is NOT equal to the counterparty- 
risk-free value, as shown in equation (26). 

Proposition 4 or equation (26) provides a theoretical explanation for the failure of full 
collateralization in the CDS market. It tells us that under full collateralization the risky value is in general 
not equal to the counterparty-risk-free value except in one of the following situations: i) the market value 
is equal to the default payment, i.e., V(u) = R(u) ; ii) firms A, B, and C have independent credit risks, i.e., 

Pij =0 and g A/iC = 0 ; or iii) p B a AC + p A a BC = d ABC . 

2.3 Numerical results 

Our goal in this subsection is to study the quantitative relationship between the CDS premium and 
the credit qualities of the counterparties and reference entity, including the default correlations and 
comrelation in a CDS contract. 

In our study, we choose a new 5-year CDS with a quarterly payment frequency. Two counterparties 
are denoted as A and B. Counterparty A buys a protection from counterparty B. All calculations arc from 
the perspective of party A. By definition, a breakeven CDS spread is a premium that makes the market value 
of a given CDS at inception zero. Assume that the reference entity C has an “A+200bps” credit quality 
throughout this subsection. The 5-year - counterparty-risk-free CDS premium is 0.027 (equals the 5-year ‘A’ 
rated CDS spread in Table 2 plus 200 basis points). 
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Since the payoffs of a CDS arc mainly determined by credit events, we need to characterize the 
evolution of the hazard rates. Here we choose the Cox-IngersoU-Ross (CIR) model. The CIR process has 
been widely used in the literature of credit risk and is given by 

dh t = ci(b - h t )dt + a-JJ^dW t (27) 

where a denotes the mean reversion speed, b denotes the long-term mean, and a denotes the volatility. 

The calibrated parameters are shown in table 4. We assume that interest rates are deterministic and 
select the regression-based Monte-Carlo simulation (see Longstaff and Schwartz (2001)) to perform risky 
valuation. 

We first assume that counterparties A, B, and reference entity C have independent default risks, i.e., 
p Alj = p AC = p BC = cp AH = Cmsc: = 0, and examine the following cases: i) B is risk-free and A is risky; and ii) 
A is risk-free and B is risky. We simulate the hazard rates using the CIR model and then determine the 
appropriate discount factors according to Proposition 3. Finally we calculate the prices via the regression- 
based Monte-Carlo method. The results are shown in Table 5 and 6. 

From table 5 and 6, we find that a credit spread of about 100 basis points maps into a CDS premium 
of about 0.4 basis points for counterparty A and about -0.7 basis points for counterparty B. The credit impact 
on the CDS premia is approximately linear. As would be expected, i) the dealer’s credit quality has a larger 
impact on CDS premia than the investor’s credit quality; ii) the higher the investor’s credit risk, the higher 
the premium that the dealer charges; iii) the higher the dealer’s credit risk, the lower the premium that the 
dealer asks. Without considering default correlations and comrelations, we find that, in general, the impact 
of counterparty risk on CDS premia is relatively small. This is in line with the empirical findings of Arora, 
Gandhi, and Longstaff (2009). 

Next, we study the sensitivity of CDS premia to changes in the joint credit quality of associated 
parties. Sensitivity analysis is a very popular way in finance to find out how the value and risk of an 
instrument/portfolio changes if risk factors change. One of the simplest and most common approaches 
involves changing one factor at a time to see what effect this produces on the output. We are going to 
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examine the impacts of the default correlations p AB , p AC , p BC , and the comrelation £ ABC separately. 


Assume that party A has an ‘A+lOObps’ credit quality and party B has an ‘A’ credit quality. The 5-year 
risky CDS premium is calculated as 0.02703. 

Assume <p AB = 0.5. The impact diagrams of the default correlations and comrelation are shown in 
Figure 2. From this graph, we can draw the following conclusions: First, the CDS premium and the default 
co rre I a t i o n s/c o m relation have a negative relation. Intuitively, a protection seller who is positively correlated 
with the reference entity (a wrong way risk) should charge a lower premium for selling credit protection. 
Next, the impacts of the default correlations and comrelation are approximately linear. Finally, the 
sensitivity slopes of the CDS premium to the default correlations and comrelation are -0.06 to p AB ; -0.09 
to p AC ; -53 to p BC ; and -14 to £ AHC ■ Slope measures the rate of change in the premium as a result of a 
change in the default dependence. For instance, a slope of -53 implies that the CDS premium would have 
to decrease by 53 basis points when a default correlation/comrelation changes from 0 to 1. 

As the absolute value of the slope increases, so does the sensitivity. The results illustrate that p BC 
has the largest effect on CDS premia. The second biggest one is £ AHC . The impacts of p AB and p AC are 
very small. In particular, the effect of the comrelation is substantial and has never been studies before. A 
natural intuition to have on CDS is that the party buying default protection should worry about the default 
correlations and comrelation. 


3 Conclusion 

This article presents a new valuation framework for pricing financial instruments subject to credit 
risk. In particular, we focus on modeling default relationships. Some well-known risky valuation models in 
the market can be viewed as special cases of this framework, when the default dependencies are ignored. 

To capture the default relationships among more than two defaultable entities, we introduce a new 
statistic: comrelation , an analogue to correlation for multiple variables, to exploit any multivariate statistical 
relationship. Our research shows that accounting for default correlations and comrelations becomes 
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important, especially under market stress. The existing valuation models in the credit derivatives market, 
which take into account only pair-wise default correlations, may underestimate credit risk and may be 
inappropriate. 

We study the sensitivity of the price of a defaultable instrument to changes in the joint credit quality 
of the parties. For instance, our analysis shows that the effect of default dependence on CDS premia from 
large to small is the correlation between the protection seller and the reference entity, the comrelation, the 
correlation between the protection buyer and the reference entity, and the correlation between the protection 
buyer and the protection seller. 

The model shows that a fully collateralized swap is risk-free, while a fully collateralized CDS is 
not equivalent to a risk-free one. Therefore, we conclude that collateralization designed to mitigate 
counterparty risk works well for financial instruments subject to bilateral credit risk, but fails for ones 
subject to multilateral credit risk. 


Appendix 


Proof of Proposition 1. Let t = T 0 . On the first cash flow payment date 7j s let V (7j ) denote the 
market value of the instrument excluding the current cash flow X l . According to Equation (10), we have 

m = E[K(T 0 ,T l ){x 1 +V(T l ))\j r t ] (Al) 


Similarly, we have 


V(T t ) = E 


K(T,,T 2 )(X 2 +V(T 2 ))n 


(A2) 


Note that K(T 0 ,T X ) is ./( -measurable. According to taking out what is known and tower 

properties of conditional expectation, we have 

V(t) = e[k(T 0 ,7j)(x, + V (7j ))|_/}]= e[k(T 0 ,T l )X l \f t ] 

i\K(T 0 ,7j )[i<XK(T\ , T 2 )X 2 \jF Ti ) + E(K(T { , T 2 )V(T 2 )| JF Ti ) 


+ E\ 


I ft 


(A3) 


: E[K(T 0 ,T l )X l \f t ]+ £[(n‘ o K(Tj , T j+1 ))x 2 \jT t ]+ £ [fn' 0 K(Tj , T j+l ))v(T 2 )\f t 
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By recursively deriving from 7 2 forward over T m , where V(T m ) = X m , we have 


m = TIA^o K ( T x T jM\^\ (A4) 

Proof of Proposition 3. Let t - T 0 , On the first payment date T, let V (7,) denote the market 
value of the CDS excluding the current cash flow . There are a total of eight (2 3 = 8) possible states 
shown in Table Al. The risky price is the discounted expectation of the payoffs and is given by 


V(t) = e\ [l 


(F(7i)X,>0 


(Pooo(^o^ 1 ) + Pioo(^o^i¥ 5 (^)+Poio(^o^i)^(^)+Piio(^o^ 1 )^ l; (r i ))(V(r i ) + Xi)|J £ 

+ \v(to + xo<o((Pooo( T o’ T 1 )+ Pm( T 0’ T i)9 A ( T i) + Po 1 o( T o’ T 1 )9a( T i) + Pm( T 0’ T i)9A B ( T i)Y V ( T l ) + X i)\ft 

+ {Pooi( T 0’ T i)+ Pm( T 0’ T iW B ( T i)+ Pou( T o’ T i)9b( T i)+ Pin( T o’ T i) ( PAB( ,r i)) R ( T s ,T l )\fty > ( t ’ T )} 

= e{[O(T 0 ,T l )(V(T l ) + X l ) + n (T 0 ,T } )R(T s , 7, )]| f t } 


(A5a) 

where 

O(T 0 ,T t ) = l (m)+Xl ),o^(7’o,7’i) + l (y(2 . i ) +Xi) <o^(r o ,ri) (A5b) 

<t> A T 0 T|) = [p A (T 0 , 7, )p B (T 0 , 7, )p c (7 0 ,T X ) + q A (T 0 , T { )p B (T 0 , T { )p c (7 0 T, )<p A (T t ) 

+ Pa T 0 >T )Ps (To >T )Pc (To >T¥a (T) + Pa <To>TV/b T 0 T)p c T 0 , T { )<p AB (7,) 

+ p c (7 0 ,7, )ct ab (ToT)( l - (P A T) - T T) + (p AB Ti)) (A5c) 

+ O-AC To, T )[p* To, T )(l - y> A T))+ ?b T 0 , T )(Ta T) - (p AB T))] 

+ °Tc To,T)[7 a To,T )(1 - Tt T)) + q A T 0 ,T)(t (T) “ <Pab T ))] 

+ Tbc To, T)(- 1 + Pa T,) - P ab (T) + Pa T ))}D(T,, T) 

K T„T,) = {p A T)T, )/? B T 0 Ti )Pc To T,) + p a (T 0 , T t )p B (7 0 ,T t ) Pc To, T )^ B T) 

+ p A T,.T)?s To - T )p c To.T )^ B T) + p A T 0 >T )? B T 0 TT c T 0 > T )<?ab T) 

+ Tc To > T )o- AB To T, )(i - ^b T) - P b T) + ^ab T)) (A5d) 

+ cr AC To, T )[p b To T, )(i - ^ B T))+ ^B To, T )(^ B (T) - p AB T))] 

+ <Tc To,T )[p A To, T 1 )(1 - P B (Ti )) ■+ Pa (To >T T) - ^ab T))] 

@ABC (T 0 ? T )(- 1 + PB (Ti ) - PAB (Ti ) + ^B T ))}DTo, T) 

QToT) = {pA ToT)P b To ,T)q c To T) + q A T 0 ,T)p B (T Q ,T)q c (T 0 , T)cp B (7,) 

+ 7a To >?> B (7 0 , 7> c To. L)^ b (7;) + <? A To> T)q B (7 0 ,7)^ c (7 0 , T)cp AB (7,) 

+ Pc (To . 7 ’)c A b To . r )(! - (’b T) - ^b T ) + Pab ( T i )) (A5e) 

— a ac To ^ T)\_p B (7 0 ,7)(l — p B (T l i) + q B T 0 > T){<p B T) — <p AB T))] 

- cr gc (7 0 ,7)[p A (7 0 ,7)(l — ^ B T )) + Pa (To , T,) - Pa* T ))] 

+ 0 ABC (7 0 ,7)(1 - P B Ti) + ^AB T) - Pb T ))}^To , n 
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Similarly, we have 


V(T t ) = E^O(T v T 2 )(X 2 +V(T 2 ))+D.(T l ,T 2 )R(T l ,T 2 )]\ J f ) (A6) 

Note that 0{T n ,T t ) is Jf T] -measurable. According to taking out what is known and tower properties 
of conditional expectation, we have 

V(t) = E{[O(T 0 ,T l )(X l +V(T l ))+n(T 0 ,T 1 )R(T 0 ,T l )]\J' t }=E[o(T 0 ,T l )X l \f t ] 

+ £[Q(r o , T x )R(T 0 , T { )\f t ] + Ep(T 0 , T x )E(fO(T l ,T 2 ){X 2 +V(T 2 ))+ 0(7;, T 2 )R(T l , T 2 )] | ) \j t } 

=s?., 4tn7J„ of 7 ', ^. ]+ z;., £ [fn7- °c r , • > 


By recursively deriving from T 2 forward over T m ? where V(T m ) - X m , we have 


(A8) 
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Table 1. Payoffs of a bilateral risky instrument 
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This table shows all possible payoffs at time t + At. In the case of V(t + At) > 0 , there are a total of four 


cases: i) Both A and B survive with probability p 00 . The instrument value is equal to the market value: 
V(t + At) . ii) A defaults but B survives with probability p l0 . The instrument value is a fraction of the market 
value: <p B {t + At)V(t + At) where tp B represents the non-default recovery rate. tp B =0 represents the one¬ 
way settlement rule, while <p B - 1 represents the two-way settlement rule, iii) A survives but B defaults with 
probability p 0l . The instrument value is a fraction of the market value: cp B (t + At)V(t + At) , where <p B 
represents the default recovery rate. 4) Both A and B default with probability p n . The instrument value is 
a fraction of the market value: cp AB (t + At)V(t + A?) , where <p AB denotes the joint recovery rate when both 
parties A and B default simultaneously. A similar logic applies to the case of V(t + At) < 0 . 


State 

o 

ll 

o' 

II 

>■? 

o 

ll 

II 

II 

>■? 

o' 

ll 

>■? 

Y a = 1,Y b =1 

Comments 

A and B survive 

A defaults, B survives 

A survives, B defaults 

A and B default 

Probability 

Poo 

P 10 

Poi 

Pn 

Payoff if VU+At) > 0 

V(t + At) 

cp B ( t + At )V (t + At) 

(p B ( t + A t)V (t + At) 

tp AB (t + At)V(t + At) 

Payoff if V(t + At ) < 0 

V(t + At) 

tp A ( t + A t)V ( t + At) 

cp A ( t + A t)V ( t + At) 

(p AB ( t + A t)V ( t + At) 


Table 2: Current/spot market data 
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This table displays the current (spot) market data used for all calculations in this paper, including the term 
structure of continuously compounded interest rates, the term structure of A-rated breakeven CDS spreads, 
and the curve of at-the-money caplet volatilities. 


Term (days) 

31 

91 

182 

365 

548 

730 

1095 

1825 

2555 

3650 

5475 


0.002 

0.002 

0.002 

0.004 

0.007 

0.010 


0.024 

0.030 

0.035 

0.040 

Interest Rate 

8 

7 

9 

3 

1 

2 

0.016 

9 

6 

5 

5 

Credit Spread 

0.004 

0.004 

0.004 

0.004 

0.004 

0.005 

0.005 

0.007 

0.007 

0.009 

0.010 

2 

2 

2 

5 

9 

2 

8 


9 

1 

6 

Caplet 

0.326 

0.331 

0.337 

0.350 

0.364 

0.377 

0.308 

0.247 

0.214 

0.167 

0.163 

Volatility 

7 


6 

9 

1 

3 


3 

1 

8 

4 
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Table 3. Impact of the credit quality of the fixed-rate payer on swap rates 


This table shows the effect of counterparty credit quality on swap rates. Assume p AB = cp AB = 0. The 1st 
data column represents the risk-free results. For the remaining columns, we assume that party R is risk-free 
and party A is risky. ‘A+lOObps’ represents a ‘ 100 basis points’ parallel shift in the A-rated CDS spreads. 


The results in the row ‘Difference from Risk-Free’ = risky swap rate - risk free swap rate. 


Credit Quality 

Party A 

~ 

A 

A+100bps 

A+200bps 

A+300bps 

Party B 


- 

- 

- 

- 

Swap Rate 

0.03433 

0.03445 

0.03459 

0.03473 

0.03485 

Difference from Risk-Free 

0 

0.012% 

0.026% 

0.04% 

0.052% 
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Table 4: Risk-neutral parameters for CIR model 


This table presents the risk-neutral parameters that are calibrated to the current market shown in Table 2. 


‘A+lOObps’ represents a ‘100 basis points’ parallel shift in the A-rated CDS spreads. 


Credit Quality 

A 

A+100bps 

A+200bps 

A+300bps 

Long-Term Mean ci 

0.035 

0.056 

0.077 

0.099 

Mean Reverting Speed b 

0.14 

0.18 

0.25 

0.36 

Volatility cr 

0.022 

0.028 

0.039 

0.056 
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Table 5: Impact of the credit quality of the protection buyer on CDS premia 

This table shows how the CDS premium increases as the credit quality of party A decreases. The 1st data 
column represents the counterparty-risk-free results. For the remaining columns, we assume that party B is 
risk-free and party A is risky. ‘A+100bps’ represents a ‘100 basis points’ parallel shift in the A-rated CDS 
spreads. The results in the row ‘Difference from Risk-Free’ = current CDS premium - counterparty-risk- 


free CDS premium. 


Credit Quality 

Party A 

~ 

A 

A+100bps 

A+200bps 

A+300bps 

Party B 

~ 

" 

" 

" 


CDS premium 

0.027 

0.02703 

0.02708 

0.02713 

0.02717 

Difference from Risk-Free 

0 

0.003% 

0.008% 

0.013% 

0.017% 
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Table 6: Impact of the credit quality of the protection seller on CDS premia 

This table shows the decrease in the CDS premium with the credit quality of party B. The 1st data column 
represents the counterparty-risk-free results. For the remaining columns, we assume that party A is risk-free 
and party B is risky. ‘A+lOObps’ represents a ‘ 100 basis points’ parallel shift in the A-rated CDS spreads. 
The results in the row ‘Difference from Risk-Free’ = current CDS premium - counterparty-risk-free CDS 


premium. 


Credit Quality 

Party A 

Party B 

A A+100bps A+200bps A+300bps 

CDS premium 

Difference from Risk-Free 

0.027 0.02695 0.02687 0.0268 0.02672 

0.00% -0.005% -0.013% -0.020% -0.028% 
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Table Al. Payoffs of a trilateral risky CDS 


This table shows all possible payoffs at time 7,. In the case of V(T [ ) + X i >0 where V(T l ) is the market 
value excluding the current cash flow X i , there are a total of eight ( 2 3 = 8) possible states: i) A, B, and C 
survive with probability p 000 . The instrument value equals the market value: V(T i ) + X r ii) A defaults, but 
B and C survive with probability p l00 . The instrument value is a fraction of the market value: 
<Pb( t i)( v ( t i) + A, ) where <p B , <p B , and cp AB are defined in Table 1. iii) A and C survive, but B defaults with 
probability p 010 . The instrument value is given by cp B (7j )(V(T I ) + X 1 ). iv) A and B survive, but C defaults 
with probability p 001 . The instrument value is the default payment: R(T lr Tp ■ v) A and B default, but C 
survives with probability p UQ . The instrument value is given by (p AB (T l ){V(T { ) + A,). vi) A and C default, 
but B survives with probability p lQl . The instrument value is a fraction of the default payment: 
<p B (T)R(T lr T l ). vii) B and C default, but A survives with probability p QU , The instrument value is given by 
cp B {T)R(T 0 ,T\) ■ viii) A, fi, and C default with probability p ul . The instrument value is given by 
cp AB (T)R(T q ,T x ). A similar logic applies to the case of V(T l ) + X 1 <0. 


Status 

Probability 

Payoff if V(T t ) + X t >0 

Payoff if V(T t ) + X t <0 

o 

II 

o' 

II 

o' 

II 

Pooo 

V(T t ) + X t 

V(T 1 ) + A, 

II 

03 

II 

o 

cH 

II 

o 

P 100 

(p^iViT^ + X,) 


o 

II 

II 

o" 

II 

P 010 

Wh (T\ )(y 17|)+a, ) 

(PMivvj + xJ 

II 

o 

II 

o 

II 

Tool 

RPo’TJ 

R(T 0 S t ) 

II 

03 

II 

II 

o 

Tno 

<p AB (T,){V(T,) + X,) 

( P,w(T l )(V(T l ) + X t ) 

II 

03 

II 

O 

II 

Tioi 

<p B (T)R(T u ,T t ) 

ip B (T)R(T 0 ,T l ) 

II 

II 

tN 

o 

II 

Ton 

<P B (T)R(T q ,T 1 ) 

<p B (T)R(T 0 ,T l ) 

II 

II 

II 

Tm 

<p ab (T)R(t 0 ,t i ) 

<Pab(T)R{t 0 ,t 1 ) 
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Figure 1. Impact of default correlation on swap rates 

The curve in this diagram represents the sensitivity of swap rates to changes in counterparty default 
correlation. Party A has an ‘A+300bps’ credit quality and party B has an ‘A’ credit quality. ‘A+300bps’ 
represents a ‘300 basis points’ parallel shift in the A-rated CDS spreads. 


Impact of Default Correlation 
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Figure 2: Impact of default correlations and comrelation on CDS premia 

Each curve in this figure illustrates how CDS premium changes as default correlations and comrelation 
move from -1 to 1. For instance, the curve ‘corBC’ represents the sensitivity of the CDS premium to 
changes in the correlation p BC when p Alj = p AC = ^ ABC = 0. 


Impact of Default Correlations and Comrelation on CDS Premia 



correlation or comrelation 


-♦— cor_AB 
— cor_AC 
-A— cor_BC 
■■— comr ABC 
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